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The ground state of an array of coupled, spin-half, antiferromagnetic ladders is studied using spin- 
wave theory, exact diagonalization (up to 36 sites) and quantum Monte Carlo techniques (up to 256 
sites). Our results clearly indicate the occurrence of a zero-temperature phase transition between 
a Neel ordered and a non-magnetic phase at a finite value of the inter-ladder coupling (a c — 0.3). 
This transition is marked by remarkable changes in the structure of the excitation spectrum. 
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Ground-state (GS) correlations in a quantum antifer- 
romagnet are closely related to the nature of low-energy 
excitations. This is deeply connected to the mechanism of 
spontaneously symmetry breaking and turns out clearly 
whenever competing GS's give rise to a quantum phase 
transition, leading in general to remarkable changes of 
the excitation spectrum at the critical point 

A simple model system experiencing such changes in 
the structure of low-energy excitations in correspondence 
of a quantum phase transition is the spin-half Heisenberg 
antiferromagnet on a two-dimensional array of coupled 
ladders: 
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where S r = are s = 1/2 operators on the 

sites r = (r x , r y ) of a L x L lattice with periodic bound- 
ary conditions; x = (1, 0), y = (0, 1), J r = J or J r = aJ 
(J > 0), depending on the parity of r y , and a is the 
inter-ladder coupling (see Fig. |l|). Such Hamiltonian in- 
terpolates between the Heisenberg model on the square 
lattice (a = 1) and a system of L/2 decoupled two-leg 
ladders (a = 0). In the square lattice limit, the GS has 
Neel long-range order, with a gapless excitation spectrum 
and a sizable value of the antiferromagnetic order param- 
eter The two-leg ladder, instead, has a finite triplet 
gap in the thermodynamic limit and no long-range anti- 
ferromagnetic order As a result a quantum critical 
point between a gapless, magnetically ordered phase and 
a non-magnetic GS of purely quantum mechanical nature 
is expected at a critical value of the inter-ladder coupling, 

Besides its intrinsic theoretical interest, this model has 
been recently studied due to the discovery of several com- 
pounds, such as SrCu2C>3 and (VO)2P20y, displaying 
clear signatures of a finite gap (~ J/2) in the excita- 
tion spectrum related to the underlying ladder structure 
HU- In addition, the coupled-ladder Hamiltonian (Q) 
has been also considered as a simplified model for the 
magnetism of the striped phase in Cu02 planes of hole- 
doped high—T c copper-oxides 

In this paper, by means of spin-wave (SW) theory, 



exact diagonalization by the Lanczos method and zero- 
temperature quantum Monte Carlo techniques, includ- 
ingvariational and Green function Monte Carlo (GFMC) 
H|| , we focus on the GS and the low-energy excitations 
of the coupled spin-ladder model (0). In particular, by 
means of a systematic size-scaling analysis, we will show 
how the structure of low-energy excitation spectrum pro- 
vides very clear indications of the changes of the GS state 
correlations occurring at a critical point, thus allowing 
us to put on firmer grounds the existence of a quantum 
phase transition at a c ~ 0.3. 
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FIG. 1. The array of coupled two-leg ladders. 

The simplest approach to the study the effects of zero- 
point fluctuations on the GS correlations of a quantum 
magnet is SW theory. This has turned out to be a very 
reliable approximation of the GS of spin-half systems 
whenever it has long-range antiferromagnetic order even 
in presence of strong quantum fluctuations 10 1. By the 
use of the standard Holstein-Primakoff representation of 
the spin operators we can compute the fluctuations over 
the classical solution at the leading order in 1/s. In con- 
trast to the Heisenberg model on the square lattice, the 
reduced translation symmetry along the y direction, im- 
plies the existence of two inequivalent branches of SW 



excitations, and the energy of such SW modes, 
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reads 
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withF, = [{li-llY+Atlilk+lk) 2 ] 1 ' 2 , 7* = (/3cosfc x + 
cos k y )/2, 6k = (1 — a)/2(l + a) sin k y , k — k + (0, tt), 
D = (1+j9)/2, and j3 = 2/(l+a). The dispersion relation 
of the two SW modes is plotted in Fig. @. Notice that 



1 



for a < 1 a finite gap develops between the optical and 
the acoustical branch of the SW dispersion along the y- 
direction due to the reduced translation symmetry. The 
acoustical branch remains always gapless at k = (0, 0) 
and k = (tt, 0) , corresponding to the Goldstone modes 
(in a reduced-zone scheme) associated with the SU(2) 
symmetry-breaking assumption. 
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FIG. 2. Linear SW dispersion relation for a = 1.0 (contin- 
uous line), a = 0.8 (short dashes), and a = 0.6 (long dashes). 

The expansion of the staggered magnetization at the 
first order in 1/s indicates the stability of the Neel order 
up to very small values of the inter-ladder coupling (a c — 
0.01, for s = 1/2). Within the same approximation the 
GS energy per site reads 
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where J — (1 + a) J/2 and the summation is over the 
fc-vectors belonging to the sector of the Brillouin zone 
with — 7r/2 < k y < 7r/2. The SW energy in the thermo- 
dynamic limit as a function of the inter-ladder coupling 
a is plotted in Fig. |^-a. For spin models the exact value 
of the GS energy can be calculated numerically for small 
clusters with the Lanczos technique and for rather larger 
sizes, in absence of frustration, with the GFMC method 
H[| . The comparison with the extrapolation to the bulk 
limit of the GFMC results (up to L < 16, see below 
for details) indicates that the latter analytical approach 
provides accurate estimates of the GS energy only in the 
regime of large inter-ladder coupling (a > 0.6) thus sug- 
gesting that the SW approach underestimates the actual 
effect of quantum fluctuations and, therefore, the critical 
value of the inter-ladder coupling. 

By means of the SW analysis it is also possible to de- 
rive a variational wave function providing a good repre- 
sentation of the GS and low-lying excited states at least 
for a — > 1, when long-range antiferromagnetic correla- 
tions are expected. This wave function, which is also 
easily computable when used for importance sampling in 
a GFMC calculation to reduce the numerical effort ||[| , 
can be obtained starting from a Neel ordered state and 
including Gaussian fluctuations by means of a Jastrow 



factor 0: 
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Here \x) is an Ising spin configuration specified by as- 
signing the value of S* for each site, Vs is the projec- 
tor onto the subspace with S^ ot = ^2 r — S, and 
Sm(x) = (— l) N i( x ) is the Marshall sign, reproducing 
exactly the phases of the GS |l2|. This depends only 
on the number of up spins on one of the two sublat- 
tices, iVf(x), so that L/V) = J2 X Sm(x)\x) represents the 
classical Neel state pfj. For the two-body Jastrow po- 
tential, the simple form, based on the consistency with 
linear SW theory in the square lattic e case [[Till, v(r) = 

(v/l 2 ) E ¥0 e~ tk - r v k , with Vk = i-V(i + r fe )/(i-r fe ), 

T/j = (cos k x + cosfc y )/2, and rj variational parameter, 
provides good variational estimates only for values of a 
very close to 1 (see Fig. ||-b), as it is expected since the 
modulation of the exchange interaction weakens the an- 
tiferromagnetic ordering and enhances spin fluctuations. 
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FIG. 3. (a): a-dependence of the exact GS energy per 
site for L — 4 (triangles), 6 (squares), 10 (pentagons), and 
12 (circles). Stars and continuous line are bulk-limit extrap- 
olations and the dashed line is the linear SW prediction, (b): 
accuracy of the variational energy of the SW-like wave func- 
tion of Eq. (^) with the isotropic Jastrow potential of Ref. |ll]] 
(empty symbols) and with its generalization to coupled lad- 
ders (full symbols). 6x6: squares; 12 x 12: circles. 

Fairly better variational results can be obtained by 
generalizing the latter wave function to the case of cou- 
pled spin ladders using the results of the SW analysis 
presented above. Indeed, the resulting wave function, 
whose explicit expression is too cumbersome to be re- 
ported here, provides very good variational estimates up 
to remarkably smaller values of a ~ 0.6 (Fig. [|-b). De- 
creasing further the inter-ladder coupling the accuracy 
of the SW variational wave function rapidly decreases. 
This clearly indicates that entering the small-a regime 
our gapless Neel ordered variational ansatz does not re- 
produce correctly the GS correlations. 
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FIG. 4. Size scaling of the GS energy per site for a = 0.1 
and a = 0.6. The dashed line is the linear fit of the data for 
L < 10 and the continuous line is the fit according to Eq. (JE]). 

Useful indications on the nature of the thermodynamic 
GS can be obtained numerically by studying the finite- 
size scaling of the GS energy, which is deeply connected 
to the nature of the excitation spectrum and therefore to 
possible thermodynamic broken symmetries [l3| ]. In fact, 
in presence of Neel long-range order, being the spectrum 
gapless and the magnon dispersion relation linear in the 
wave vector k, the leading finite-size correction to the GS 
energy per site, e Q (L) = E Q (L)/L 2 , is 0(L' 3 ) gj]- In- 
stead, in presence of a finite correlation length, a finite 
gap in the spin excitation spectrum and an exponential 
asymptotic dependence is expected. Here we have as- 
sumed a size dependence of the form 



e (L) = e (oo) + c exp(-L/L )/L 2 



(5) 



which has been employed in Ref. jl5[] for the two-leg lad- 
der. As it is shown in Fig. ^, the size-scaling law predicted 
for a long-range ordered GS is fulfilled for values of a 
close enough to 1 while for small value of a a clear devi- 
ation from this behavior is observed and the exponential 
law (||) is instead satisfied. This provides another nu- 
merical evidence of the melting of the antiferromagnetic 
long-range order due to the transition to the ladder-like 
regime. 

The occurrence of a quantum phase transition to a non- 
magnetic GS is clearly confirmed by the study of the spin 
gap shown in Fig. |^. This physical quantity can be mea- 
sured straightforwardly with GFMC by performing two 



different simulations in the S? 



and Sf. 



1 sub- 



spaces. In contrast to finite-temperature algorithms, the 
calculation of the spin gap within GFMC does not in- 
volve any fitting procedure of low-temperature data and 
it is exact within statistical error. As shown in the left 
panel of Fig. [j| the triplet gap, A = E\ — Eq, is a de- 
creasing function of the inter-ladder coupling even if, on 
finite-sizes, this quantity is always non-zero for any value 
of a. However, the size scaling of the spin gap, shown 
in the right panel of the same figure, indicates a clear 



deviation from the size dependence expected in presence 
of long-range Neel order E3] , 



A L =a/L 2 + b/L 3 



(0) 



and the opening of a finite gap in the thermodynamic 
excitation spectrum for a < 0.35. In this regime the 
extrapolation to the bulk limit of the finite-size data can 
be done using a law of the type |16| : 



A L = A + a/L 2 + b/L 4 



(7) 




FIG. 5. Left panel: a-dependence of the triplet gap for 
various lattice sizes. Stars are bulk-limit extrapolations ac- 
cording to Eq. (M) and the dashed line is a fit according to 
A oc (a c — a) '. Right panel: size scaling of the triplet gap 
for different values of a. Continuous and dashed lines are fits 
according to Eqs. (j?|) and respectively. 

An estimate of the critical value of the inter-ladder 
coupling can be obtained by fitting the extrapolated val- 
ues of the gap with the scaling law A = (a c — a) v , 
with v ~ 0.69, predicted for a quantum phase transi- 
tion in (2+1) dimensions This procedure gives the 
value a c ~ 0.32 ± 0.03, in agreement with previous nu- 
merical estimates obtained with finite-temperature algo- 
rithms (|;@] an d t ue mean-field predictions of Ref. fl7f| 
but in contrast with the conclusions of bond-mean-field 
theory J(| indicating the vanishing of the gap for infinites- 
imal values a. In contrast to frustrated systems like the 
J\ — Ji model, where the SW theory provides an accurate 
prediction of the transition to a non-magnetic phase due 
to competing interactions Jl0| , ^8 19 1, in this case the SW 
result (a c ~ 0.01) grossly underestimate the exact one. 
This can be ascribed to the fact that in this model there 
are no competing GS's at the classical level, the Neel 
state being stable up to a = 0, and the transition has 
therefore a pure quantum origin. Within the 1/s expan- 
sion, instead, the SW velocity remains finite up to a = 
and the reduction of the staggered magnetization is only 
due to the crossover to the one-dimensional regime in 
which antiferromagnetic long-range order is unstable. 
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FIG. 6. Lower panels: rigid rotator anomaly vs S for 
a = 0.2, and a = 0.8: L — 6 (triangles), 8 (squares), 10 
(circles). Upper panel: size scaling of Sx~ 1 (l)/(L/2 - l)L 
for (from the top) a=0.2, 0.25, 0.3, 0.4, and 0.8. Lines are 
weighted quadratic fits. 

A major fingerprint of the dramatic changes in the 
structure of the excitation spectrum occurring in corre- 
spondence of the above quantum phase transition can 
be found in the finite-size behavior of the so called rigid 
rotator anomaly |po|, Sx~ 1 (S) = l/2\s ~ 1/2xl, with 
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where Eg is the energy of the lowest excitation with spin 
S. In fact, in presence of long-range antiferromagnetic or- 
der, the low-lying excited states of spin S are predicted 
to behave as the spectrum of a free quantum rotator, 
E s -E oc S{S + 1)/L' 2 , as long as S < L 2 f|. In con- 
trast, the expected behavior for a spin ladder is Eg— En oc 
S, as it is easy to understand in a spin-liquid Resonat- 
ing Valence Bond (RVB) picture . As a result, 
in the gapless phase the rigid rotator anomaly 5\ 1 has 
to vanish identically in the thermodynamic limit while 
in the gapped regime 5x~ 1 (S) oc L(L/(S + 1) — 1) di- 
verges linearly with the volume. These features of the 
excitation spectrum clearly discriminate the two differ- 
ent zero-temperature phases of the present model, as it 
is evident from Fig. || displaying the very different be- 
havior of Sx~ 1 below and above the critical point. 

In summary, we have investigated the ground-state 
properties of an array of coupled spin-half ladders using 
spin-wave theory and a numerical analysis of the finite- 
size low-energy excitation spectrum. While the spin- wave 
theory turns out to be reliable only in the regime of weak 
inter-ladder coupling, our numerical results provide ro- 
bust indications of the opening of a finite spin gap in 
the thermodynamic limit for a < 0.3, corresponding to a 
quantum phase transition between a gapless Neel ordered 
phase and a spin-liquid RVB ground state. 
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